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Abstract. The steady state velocity fluctuations of a movable piston located on the 
top of a vibrated granular gas are studied by means of molecular dynamics simulations. 
From the second moment of the distribution, a temperature parameter for the piston is 
defined and compared with the granular temperature of the gas located just below it. 
Then the two temperature parameters refer to two interacting macroscopic systems. 
The equipartition of energy valid in usual molecular systems is strongly violated, 
and the temperature of the piston can be larger or smaller than the one of the gas, 
depending of the parameters defining the system at the particle level. The simulation 
results for the ratio of temperatures are in agreement with some theoretical predictions 
from kinetic theory assuming the validity of a hydrodynamic description in the limit 
of weak inelasticity of the gas. 
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One of the characteristic features exhibited by granular fluids is the violation of the 
energy equipartition theorem. The granular temperatures of each of the components of 
a mixture, defined from the average kinetic energies of the particles, are not the same. 
This was pointed out long ago [1] , and has attracted a lot of attention in the last decade 
or so. It is important to stress that the equality of the temperatures of the components 
of a mixture of molecular fluids is not restricted to equilibrium thermal systems, but 
also applies to out-of-equilibrium situations, at least to those that can be described 
by means of usual hydrodynamics as derived from kinetic theory or non-equilibrium 
statistical mechanics. 

For the homogeneous cooling state (HCS) of a binary granular mixture modeled 
as an ensemble of smooth inelastic hard spheres, the ratio of temperatures of the 
two components has been obtained from an approximated solution of the kinetic 
Enskog equations [2]. This result is in good agreement with Molecular Dynamics 
(MD) simulation data [3]. The HCS is an idealized state in which the system cools 
homogeneously and monotonically. In real experiments, to maintain a granular system 
fluidized, external energy must be continuously supplied, and this is achieved by means 
of vibrating or rotating walls, by shearing, by shaking, or by other external fields. As 
a consequence, the systems become inhomogeneous. The non-equipartition has been 
confirmed by MD simulations of simple shear flows [3], and of vibrated gases [5l [6]. 
There are also some experimental evidences of the coexistence of different component 
temperatures in strongly vibrated granular mixtures [71 |8]. Homogeneously driven 
granular fluids submitted to stochastic forces have also been considered. The lack of 
energy equipartition in these systems has been analyzed [9] , but the relationship between 
this kind of ideal driving and actual experiments is uncertain. 

A peculiarity of the above studies is that they are concerned with quantities 
going beyond the macroscopic description, in the sense that the relevant macroscopic 
local property is the global temperature (or energy density) of the mixture, and not 
that of each of the components. Moreover, the latter are rather difficult to measure 
experimentally. Here, a possible experimental set-up in which the breakdown of energy 
equipartition is macroscopically observed is described and studied by means of MD 
and kinetic theory. The experiment is specifically designed to be reproducible in the 
laboratory, and might have some practical advantages over previous related experiments 
[To] , although it must be stressed that the study reported here is restricted to low density 
gases and that no friction between the particles, of the particles with the piston, or of 
the latter with the walls of the container is taken into account. 

Consider a granular system composed of inelastic hard spheres {d = 3) or disks 
{d = 2) of mass m and diameter a. There is an external gravitational field acting on 
the system, so that each particle is submitted to a force —rngQCz, where go is a positive 
constant and the unit vector in the direction of the positive z axis. The system is 
kept fiuidized by vibrating the lower wall of the container in a sawtooth way, with very 
high frequency, negligible amplitude, and velocity Vw = vw^z [HI [12]. The latter will 
be always taken large enough as to keep the density of the system small everywhere. 
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Figure 1. Sketch of the system considered in this work. 



On the top of the system there is a movable hd or piston of mass M, as illustrated in 
Fig. [H The piston can only move in the 2;-direction, remaining always perpendicular to 
it. It is assumed that there is no friction between the piston and the lateral walls of the 
vessel containing the granular gas. 

Collisions between particles are inelastic and characterized by a constant coefficient 
of normal restitution a, so that when a pair of particles i and j collides, their relative 
velocity Vij changes instantaneously according with 

Vij v[j = Vij - (1 + a)a ■ Vija, (1) 

while the total momentum is conserved. In the above expression, is a unit vector 
directed from the center of particle j to the center of particle i at contact. Collisions 
of particles with the piston are also modeled as smooth and inelastic, being ap the 
coefficient of restitution in this case. Taking into account that the piston can only move 
in the z direction, the change in the relative velocity in a collision is restricted to the 
2;-component, and it is given by 

v,-V,^v',- VI = v,-V,-{l + ap){v, - K), (2) 

where Vz is the velocity of the piston. This relation and the total momentum 
conservation define the collisions with the piston. 

Next, some results from MD simulations of two-dimensional systems will be 
presented. To avoid undesired effects, periodic boundary conditions are used in the 
direction perpendicular to the external field. It is observed that, at least when the 
width W of the system is not too large, a steady state with no macroscopic fiow of 
mass, and gradients in the fiuid only in the direction of the external field, i.e. the z axis, 
is reached after some transitory period. This state is independent from the value of W , 
as long as it is kept small enough. The shape of the hydrodynamic profiles in the bulk 
will be analyzed elsewhere IHj, while here the emphasis will be put on the relationship 
between the temperature parameters of the piston, Tp, and of the gas in its vicinity, 
Tq- Both temperatures are defined from the second moment of the respective velocity 
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distributions, namely M{V^) = Tp, m(t>^) = 2Tg, with the angular brackets denoting 
average. In the results to be reported in the following, units defined by m = 1, a = 1, 
and go = I will be used. 

The measured probability distribution for the velocity of the piston is very well 
fitted by a Gaussian in all the cases considered. A typical example, corresponding to 
vw = 11, a = 0.9, ap = 1, M = A9, N = 420, and W = 70, is given in Fig. H For this 
case, the value Tp ~ 58.67 is found. It is seen that the Gaussian fit is accurate at least 
up to values of the probability density of the order of 10~^. The measurement of Tq is 
somewhat more problematic, since the 'top region' of the granular gas next to the piston 
has to be identified. Here, the layer of the gas defined by (Z) — 2az < z < {Z) + 2az, 
where az is the dispersion of the values of the position of the piston Z, has been 
used. Then, the average over the velocities to compute Tq has been carried out over 
the particles inside this region at each time. This criterium may appear as rather 
arbitrary, but it has been checked that other sensible choices, such as extrapolating the 
bulk temperature profile to (Z), lead to values of the temperature with discrepancies 
of the order of 5%. In this way, it is found that Tq — 37.69 , implying Tp/Tc — 1.56. 
This represents a strong violation of the energy equipartition, specially taking into 
account that the collisions between the grains and the piston have been modeled as 
elastic (ap = 1). Moreover, it must be stressed that the temperature parameter of 
the piston is larger than the temperature of the gas transmitting to it the effects of 
the lower vibrating wall. It can be wondered whether this effect is related with some 
strong anisotropy of the velocity distribution of the gas. Let us define Tqx = m < > 
and Tqz = m < v"^ >, so that Tq = (Tgx + Tgz)/2. Then, for the particular case we 
are discussing, it is found that Tgz/Tgx — 1.096. Therefore, although there is some 
anisotropy as usual in vibrated granular systems, it is very weak as compared with the 
excess value of Tp/Tc- More precisely, it is Tp/Tcz — 1-49. The same trend is observed 
in all the other cases investigated, being the above the one in which the anisotropy of 
the second moment of the velocity tensor is stronger. 

In Fig. [3], the temperature ratio Tp/Tc is plotted as a function of the coefficient 
of restitution of the gas a for three different values of the pair (M, ap), illustrating 
the influence of each of these parameters on the temperature ratio. Quite interestingly, 
the violation of equipartition (deviation of the temperature ratio from unity) is not a 
monotonic function of a, but the temperature ratio exhibits a minimum in the small 
inelasticity region. The result with a = 1, ap = 1, has been obtained with a non- 
vibrated system at thermal equilibrium. As expected, equipartition is recovered in the 
limit a = ap = 1. The reason to restrict the analysis to the inelasticity region a > 0.9 
is that, for smaller values, even the hydrodynamic description of the bulk of the gas, as 
provided by the Navier-Stokes equations, breaks down [13]. 

An expression for the temperature ratio can be derived as follows. Given the low 
density of the gas, a pair of coupled Boltzmann and Boltzmann-Lorentz kinetic equations 
are used as the starting point to describe both the gas and the piston. The subsequent 
calculations are based on the following two assumptions: (a) the macroscopic state of the 
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Figure 2. Steady velocity distribution of the piston P{Vz) in both logarithmic and 
normal scales. The symbols are from MD simulations while the solid line is a Gaussian 
fitting. The values of the parameters defining the system and the dimensionless unit 
for the velocity Vz are given in the main text. 




1.6 



o o M=49 a =1.0 
m-m M=14 a=1.0 
A A M=49 a =0.9 



1.4 



1.2 



0.9 



0.92 0.94 



0.96 



0.98 



Figure 3. Ratio of the temperature of the piston, Tp, to the temperature of the gas 
next to it, Tq, versus the coefficient of restitution of the gas, a. 



piston can be specified in terms of tlie liydrodynamic fields of tlie gas in its vicinity and 
its position probability distribution, and (b) the associated normal solution of the kinetic 
equations can be generated by the Chapman-Enskog procedure. With these hypothesis, 
it is obtained that the zeroth order in the gradients cooling rates of the piston, (p^ and 
of the top region of the gas, (q \ must be the same [TH [13]. These rates are functionals 
of the one-particle distributions for the top region of the fluid and the piston. Given that 
these distributions are not known, Maxwellians can be used as an estimation. Therefore, 
no anisotropy effects are introduced in the theory. A generalization of the derivation in 




Figure 4. Ratio </> = mTp / MTq of the tlicrmal velocities for the piston and the gas 
just below it. The Hnes are from Eq. ([6]) and the symbols MD simulation results for 
hard disks. 



|13] . with the pecuharity that now all the the collision vectors have the same direction, 
leads to pjj 

0-^-^(1 + ^)"^ (3) 

ucVtct'^-^ ~ r (d/2) d 



^(0) /2 (d-l)/2 



where hq is the gas density in the top region, vt = {2TG/mY^'^ , = mTp/MTc, and 
_ m{l + ap) 
^ = 2(m + M)- 

Note that all the dependence on the mass ratio and the inelasticity ap occurs through 
h. Requiring the cooling rates to be the same, an equation determining is obtained, 

^ ^ V <P J 2V2r{d/2)Wd ^ ' 

The elastic collisions limit {a = ap = 1) is given hy cj) = h/(l — h) = m/M, that is the 
equipartition result. If the particle collisions are elastic but not those with the piston 
(a = 1, ap < 1), it is = m(l + ap)/[2M + m(l — ap)]. This result is the same as the 
one derived in ref. |l5j for an inelastic intruder in an equilibrium elastic gas. 

The prediction made by Eq. ([6]) is compared with some MD simulation results in 
Fig. m where is plotted as a function of h. In the simulations, the values of the 
latter quantity have been modified by changing both M (in the range 10-140) and ap 
(between 0.6 and 1). The other parameters {W, vw, and A^) are kept the same as in 
Fig. [21 The symbols are the simulation results while the lines are from Eq. ([6]), as 
indicated. Data for two values of a, 0.98 and 0.9, are reported. For the most elastic 
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Figure 5. Marginal probability distribution of the z-component of the scaled velocity 
of the gas, = v^/vt, in the top region. The values of all the parameters are the 
same as in Fig. [5] 

gas, the agreement is quite good, for all the considered values of the other parameters. 
Nevertheless, when the inelasticity of the gas is slightly increased, strong deviations from 
the theoretical prediction show up. These discrepancies are not only quantitative, but 
also qualitative. The circled symbols on the left lower corner of the figure correspond 
all them to M = 49, differing in the value of ap. It is observed that when ap (and 
h) increases, decreases, in contradiction with the solution of Eq. ([6]). To identify 
the cause of the failure of the theory, the shape of the marginal distribution of the 
component of the velocity of the grains perpendicular to the piston in the top region, 
ip{vz), has been measured. While for a = 0.98 it is always well fitted by a Gaussian 
within the numerical error of the simulations, for a = 0.9, exponential tails are clearly 
identified for positive velocities already in the thermal range, as it is seen in Fig. [5l 
These tails are relevant for the lowest velocity moments of the gas. On the other 
hand, it was already mentioned above (see Fig. [2]) that the velocity distribution of the 
piston remains Gaussian. A theoretical analysis of the existence of this exponential tail 
probably requires going beyond the Chapman- Enskog method, the restriction to normal 
solutions of the Boltzmann equation and, therefore, of the hydrodynamic description. 

The above results seem to provide additional evidence of the difficulties when trying 
to give a macroscopic meaning to the granular temperature, at least as defined in kinetic 
theory and the usual hydrodynamic theories. The knowledge of the hydrodynamic 
profiles of the confined gas and the characteristics of the interactions between the gas 
particles and the piston is not enough to determine the "temperature" of the latter. 
More detailed information about the velocity distribution of the gas in its vicinity is 
required. Moreover, it is surprising that sometimes the "temperature" of the body to 
which energy is supplied (the piston) be larger than the temperature of the system 
supplying the energy (the neighbor gas). The results also show that this effect is 
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not associated to the anisotropy of the gas. Finally, let us stress the differences of 
this system as compared with a gas mixture, even in the tracer (one-intruder) limit. 
Here, two macroscopic systems spatially differentiated are interacting at contact, their 
respective states being strongly interdependent. As a consequence, it seems clear that 
both temperature parameters have a macroscopic meaning. 
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